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ABSTRACT 

Recently discovered exoplanets on close-in orbits should have surface temper- 
atures of 100's to 1000's of K. They are likely tidally locked and synchronously 
rotating around their parent stars and, if an atmosphere is absent, have sur- 
face temperature contrasts of many 100's to 1000's K between permanent day 
and night sides. We investigated the effect of elevated surface temperature and 
strong surface temperature contrasts for Earth-mass planets on the (i) pattern 
of mantle convection, (ii) tectonic regime, and (iii) rate and distribution of par- 
tial melting, using numerical simulations of mantle convection with a composite 
viscous/pseudo-plastic rheology. Our simulations indicate that, if a close-in rocky 
exoplanet lacks an atmosphere to redistribute heat, a > 400 K surface tempera- 
ture contrast can maintain an asymmetric degree 1 pattern of mantle convection 
in which the surface of the planet moves preferentially toward subduction zones 
on the cold night side. The planetary surface features a hemispheric dichotomy, 
with plate-like tectonics on the night side and a continuously evolving mobile 
lid day side with diffuse surface deformation and vigorous volcanism. If volcanic 
outgassing establishes an atmosphere and redistributes heat, plate tectonics is 
globally replaced by diffuse surface deformation and volcanism accelerates and 
becomes distributed more uniformly across the planetary surface. 

Subject headings: planetary systems — planets and satellites: interiors — planets and satel- 
lites: surfaces — planets and satellites: tectonics. 

1. Introduction 



Recent discoveries of ex oplanets with 
Earth-like mass and radius ( iMayor et al. 
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2009; iHoward et al. I 120101 : iBorucki et al. 
20 111 ) have intensified debate on how such 
planets compare to Earth in other respects. 
Due to observational bias, many of the 
discovered exoplanets inhabit short-period, 
close-in orbits and have effective tempera- 
tures exceeding many hundreds of Kelvin. 
For such planets, a likely outcome of dynam- 
ical evolution is tidal locking to their par- 
ent stars and, most proba bly, capture into 



a 1:1 spin-orbit resonance (ICorreia fc Laskar 
20101 ). Synchronous rotation causes asym- 
metric insolation and, in the absence of 
a substantial atmosphere, a strong (100's 
to 1000's K) temperature contrast between 
these planets' permanent day and night 
side s, as has been estimated for CoRoT- 



7b ( iLeger et 



al. 



2009 



f) and is p lausible for 



Kepler- 10b (IBatalha et al. 11201 If ). 



Characterization of exoplanet surfaces is 
challenging and numerical simulations can 
help determine possible scenarios. Subsolidus 
convection is likely within the silicate mantles 
of rocky exoplanets but its vigor and surface 
expression depend on mantle temperature, 
composition, and rheology. Within the So- 
lar System, Earth is the only planet that cur- 
rently exhibits plate tectonics which, on a ge- 
ologic time scale, regulates volatile species in 
Earth's atmosphere via volcanic outgassing, 
silicate weather ing, and subduction of precip- 



itated carbon (IWalker et al. I Il98lh . Venus 



is thought to experience infrequent global- 
scale resurfacing events, possibly the re- 
sult of mantle-wide episodic overturn and 
associated with the formation of a dense 



greenhouse atmosphere, e.g.. ISolomon et al 



(119991 ) . Mars has experienced recent sub- 
stantial volcanism although there is no evi- 
dence for recent crustal mobility, suggesting 
a stagnant lid regime with convective activity 
in the mantle interior below a static single- 



plate lithosphere, e.g. ISpohn et al. I ( 1200 ll ). 



Elevated surface temperatu res can affect a 



plane t's interior dynamics ( jLenardic et al. 



20081 ) and its surface properties, as exempli- 
fied by insolation-driven va riations of Mer- 



cury' s lithospheric strength ( [Williams et al. 



20111 ). Although there are no examples in our 
Solar System, a hemispheric contrast in sur- 
face temperature may influence the interior 
dynamics and surface expression of exoplan- 
ets. 

Here, we investigate the effect of elevated 
surface temperature and strong surface tem- 
perature contrasts for Earth-mass planets on 
(i) patterns of mantle convection, (ii) tectonic 
regimes, and (iii) the rate and distribution 
of partial melting (volcanism). For this pur- 
pose, we conducted numerical simulations of 
mantle convection with imposed surface tem- 
peratures. 



2. Description of Numerical Model 
and Parameters 

We investigated planetary mantle con- 
vection through numerical simulations of an 
incompressible fluid at infinite Prandtl num- 
ber using the classical Boussinesq formu- 
lation. The associated conservation equa- 
tions for mass, momentum, and energy are 
solved numerically using the finite element 
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package CitComS-3.1.1 flZhong et al. Il2000h . 
We investigated effective Rayleigh numbers 
Ra eff = p g a ATh 3 /(r] eff K ) in the range 
~10 5 — 10 7 , where p , g , a , and k are the re- 
spective reference values for density, gravita- 
tional acceleration, thermal expansivity, and 
thermal diffusivity, r] e ff is the time- averaged 
mantle viscosity, AT is the temperature con- 
trast across the mantle, and h is the man- 
tle thickness. Bottom and internal heating 
both contribute to the heat budget. Inter- 
nal radiogenic heat production is uniformly 
distributed across the mantle domain and 
constant in time in our models. The non- 
dimensional internal heat genera tion rate, de- 
fined as 7 = h 2 Ho/(KCpAT) ( jGlatzmaier 



the lateral and radial directions, respectively, 
with gradual mesh refinement towards the 
top and bottom boundaries where tempera- 
ture variations are generally greatest. 

To permit plate-like behavior of the 
surface boundary layer, we adopt a com- 
posite viscous /pseudo- plastic r h eology in 



our models, following iTackley I ( 12000al Jbl). 



Temperature-dependent viscosity rj v is de- 
scribed by an Arrhenius-type law: 



Vv(T') 



?7oexp 



23.03 



1 



1 + T> 



(i; 



19881 ). is set to 11; this corresponds to a di- 
mensional heating rate of H Q = 4 x 10~ 12 
W kg -1 , clos e to the present-day chondrit ic 
heating rate (ITurcotte and Schubert II2002I ). 



Our calculations ignore dissipation of 
tidal forces that can contribute to mantle 
heating and thermal runaway, depending on 
orbital eccentricity and type of resonance 
( iBehounkova et al. 11201 ll ). However, without 
external perturbations, tidal dissipation ren- 
ders orbital eccentricity insignificant on Myr- 
Gyr time scales and tidal effects can be ne- 
glected. 

For the mantle domain, we employ an 
annulus (bi-section of a 3-D spherical do- 
main) that is aligned with the ecliptic plane. 
All domain boundaries are free-slip and two 
side boundaries are imposed at the antistel- 
lar point. We adopt Earth's outer radius (r s 
= 6371 km) and mantle thickness (h = 2891 
km), and allocate 257 x 65 nodal points in 



where the reference viscosity is r)o = 
5 x 10 20 Pa s, and T' is the dimension- 
less mantle potential temperature that re- 
lates to the dimensional temperature, T, as 
V = (T - T S )/AT, with AT = 2400 K 
the potential temperature contrast across the 
mantle and T$ the surface potential tem- 
perature. Viscosity changes by 5 orders of 
magnitude over the considered temperature 
range (T' G [0, 1]) and generates a litho- 
sphere over a weaker mantle. Pseudo-plastic 
yielding concentrates strain and allows for 
lithospheric break-up in confined regions that 
mimi c subduction zones and spreading cen - 



ters flTacklev I l2000al : iRichards et al. ||200jJ ). 



In regions where the model stress exceeds 
an assigned yield stress, a y , a yield viscos- 
ity is calculated as rj y = a y /2e'u, where 
e'n is the second invariant of the strain 
rate tensor. The composite rheology is de- 
scribed as rj = min (r) v (T), r) y ). We did 
not consider strain - or strain-rate-weakening 
(IChristensen 1119841 ). or time-dependent dam- 
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age rheology (jBercovici I Il996l ). which could 
further enhance plate-like behavior but would 
significantly complicate our analysis. 

We assign distinct temperature condi- 
tions for 3 contrasting cases. In a first set 
("cold" or C- models) we apply a uniform, 
time-constant surface temperature T s = 273 
K, comparable to Earth. A second set ("hot" 
or H-models) mimics close-in planets (i.e., or- 
bital distance of a = 0.13 AU around a so- 
lar mass star) with efficient heat redistribu- 
tion and a uniform surface temperature of Ts 
= 759 K. A third set ("asymmetric" or A- 
models) considers planets on a similar close- 
in orbit that lack heat redistribution. For 
these models, a high substellar temperature 
of T subst = 1073 K decreases sinusoidally to 
the terminus and is kept constant on the night 
side at Ts = 273 K. For all models, the core- 
mantle boundary (CMB) potential tempera- 
ture is uniform and constant at Tqmb = 2673 
K. Each simulation is run for several billion 
years (Gyr) of model time and we exclude 
the first 1 Gyr of initial transients to focus 
on statistically steady-state behavior. 

To quantify tectonic regimes, we make 
use of two pr eviously-defined diagnostics 
(ITackley I l2000al ) . First, to quantify the lo- 
calization of surface strain rates, we define 
"plateness" as P = 1 — (/so/0. 6), where / 80 is 
the area fraction that encompasses 80% of the 
total surface strain rate. P = corresponds 
to strain localization for isoviscous convec- 
tion. Second, lid mobility M is defined as the 
ratio of the root-mean-square (RMS) surface 
flow velocity relative to the RMS velocity of 



Models with plate-like behavior are charac- 
terized by M f« 1 — 1.5 and for stagnant 
lid convection M ~ 0. To allow faster cal- 
culations, we made use of the symmetry of 
the problem and determined tectonic regimes 
(Section for models of 180° opening angle, 
with side boundaries at the sub-stellar and 
anti-stellar points. For 4 representative mod- 
els the time-averaged P and M values differ 
by only 6% and 9%, respectively, when com- 
paring models of 180° and 360°. 

Pressure-re lease partial me lt ing is calcu- 



lated following iRaddick et al. I (120021 ) . We 



the entire mantle domain, M 



v 



srfc 

RMS I U RMS 



/V 



whole 



only consider melting in regions where con- 
vective flow is upward and where real (po- 
tential + adiabatic) temperatures exceed the 
mantle solidus temperature (T r > T so i). For 
T so i we use a parameterization for dry mantle 
peridotite with a 1 bar value of T® ol = 1373 
K and dT so i/dz = 3.3 K km -1 , in rough 
agreement with experimental measurements 
(e.g., Hirschmann 2000). In melting regions, 
the local melting rate, q m , is calculated as 
q m (x,t) = df/dt = (df /dz)u z (x,t), where 
/ is the degree of melting and u z = dz/dt 
is the upward convective flow velocity. We 
use a constant value df/dz = 0.18% km -1 
for the adiabatic m elt production per kilo - 
meter of upwelling (IPhipps Morgan I l200ll ). 
The local melting rate, q m , is subsequently 
integrated over each melt column to give the 
surficial melt distribution. We limit melt pro- 
duction to a maximum depth of 50 km, which 
results in a time-averaged melt production 
for a nominal Earth model (C150) consistent 
with Earth's pr esent-day melt production o f 
~ 20 km 3 yr _1 ^McKenzie and Bickle Ifl988l ). 
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Our choice of maximum depth affects the to- 
tal melt production but has a small influ- 
ence on the comparison between models or 
on surficial distributions of melt. Because 
side boundaries promote vertical flow and un- 
physical concentration of melt at the substel- 
lar point, we calculate melting from models 
with a 360° opening angle. 



3. Tectonic Regimes 

C-models show a variety of convection 
regimes for a progressive increase of the yield 
stress values a y . At o y < 50 MPa, con- 
tinuous yielding prevents the formation of 
stiff surface plates. Instead, a mobile lid 
style of convection, characterized by diffuse 
surface deformation, occurs (Figure [Jji, red 
curve). At o y ~ 150 MPa, surface defor- 
mation is approximately plate-like with stiff 
surface plates separated by narrow regions of 
concentrated deformation. Surface plates ex- 
hibit approximately piecewise constant veloc- 
ities (Figured^, green curve) and converge to- 
ward subduction-like downwellings while di- 
verging at localized spreading centers (Figure 
dU) . A further increase of the yield stress pro- 



duces a more time- dependent solution with 
alternating periods of surface mobility and 
stagnation. For sufficiently high yield stress 
(800 MPa), surface mobility diminishes and 
stagnant lid convection persists (Figure [Jji, 
blue curve). To further demonstrate tectonic 
regimes, we show plateness P and mobility 
M in Figure [2J For the C-models (Figured), 
plateness increases to P ~ 0.75 with increas- 
ing yield stress at high mobility (M ~ I — 1.5) 
until the lid mobility rapidly decreases (M ~ 
0) around a y ~ 250 MPa (and plateness be- 
comes irrelevant). 

For the H-models, coherent lithospheric 
plates do not develop because low viscosi- 
ties near the warm (Tg = 759 K) surface 
promote viscous deformation and suppress 
pseudo-plastic yielding, consis t ent w ith the 



predictions of iLenardic et al. I ( 120081 ). Sur- 



face velocities are greater and have a more dif- 
fuse distribution than in the C-models (Fig- 
ure [Tb). Due to increased lithospheric mo- 
bility (e.g. model H150, Figure [Tfe) plate- 
ness is consistently lower for the H-models 
(P < 0.5, Figure [2b) compared to the C- 
models (P ~ 0.75, As a y increases, stag- 
nant lid prevails but plate-like behavior is not 
observed (Figure [2b). 
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Fig. 1. — Top row: Surface velocities for (a) cold symmetric C- models, (b) hot symmetric 
H- models, (c) hot asymmetric A- models (see Section 2 for description). For all model types, 
3 cases are shown with different yield stresses a y of 50 MPa (red curve), 150 MPa (green 
curve), and 800 MPa (blue curve). Middle row (d-f): Snapshots of dimensionless potential 
temperature and convective flow velocity for models with a y = 150 MPa that correspond 
to the green curves in the top row. Bottom row (g-i): Time-averaged values of the same 
quantities and models as in the middle row. 



For the A-models, a marked (AT$ = 800 
K) temperature contrast causes global-scale 
asymmetries in tectonic regime. Near the 
substellar point, high surface temperatures 
induce convective upwelling and diffuse sur- 
face velocities (Figure [lb and [If). Away from 
the substellar point, plate-like behavior oc- 
curs and surface plates preferentially move 



toward the antistellar point (Figure [If). A 
global scale dichotomy is reflected in the tec- 
tonic regimes diagnosed separately for the 
day and night sides (Figure [2b). The hot day 
side is consistently characterized by mobile 
lid convection with diffuse deformation, simi- 
lar to the H-models, while the cold night side 
exhibits plate-like behavior for a large param- 
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eter space, similar to the C-models (Figure 
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Fig. 2.— Diagnostics of tectonic regime 
(plateness, P, and mobility, M) for mod- 
els with different Rayleigh number and yield 
stress, (a) "Cold symmetric" C-models, (b) 
"Hot symmetric" H- models, and (c) "Hot 
asymmetric" A-models, for which both day 
(right) and night (left) sides are shown, (d) 
All models in P-M domain space:C-models 
(dark purple diamonds), H- models (light pur- 
ple stars), night side of A-models (dark blue 
left semi-circles), and day side of A-models 
(light blue right semi-circles). Color coding in 
(a-c) is for P- values if M>0.75 and for M- 
values if M<0.75, as indicated by the color 
bars above frame (d). 



4. Patterns of Mantle Convection 

Flow patterns in the A-models are char- 
acterized by convective upwelling near the 
hot sub-stellar point, near-surface flow from 
the hot day side to the cold night side where 
most downwellings occur, and a deep mantle 
return flow toward the day side (Figured):). 
An asymmetric degree 1 pattern of mantle 
flow persists, with convection cells that oc- 
cupy the entire half-mantle domain (Figure 
[TJ). For the C- and H- models with uniform 
surface temperature, convective downwellings 
are more randomly distributed across the do- 
main and this results in less persistent con- 
vective flow (Figure [Tg and[Th). As a diag- 
nostic for persistent mantle flow, we use the 
RMS value of the time-averaged flow velocity 
normalized by the time average of the RMS 
flow velocities, or (3 = (v)rms/vrms, where 
overlines indicate time-averaged values. The 
persistent flow in model A150 is reflected by 
relatively high time-averaged convective flow 
velocities ((3 = 0.53), compared to models 
with a uniform surface temperature, C150 ((3 
= 0.41) and H150 (/9 = 0.43). 

To estimate the surface temperature con- 
trast ATs that is required for maintaining 
asymmetric mantle flow, we compare power 
spectra of time-integrated lateral convective 
flow velocities for models with different ATs 
(Figure [3]). For surface temperature con- 
trasts >400 K, a dominant degree 1 signal 
reflects persistent convection cells with up- 
welling at the substellar point and down- 
welling at the antistellar point. A systematic 
degree 1 signal is not discernible for models 
with ATs<400 K and this reflects more ran- 
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domly oriented convection cells with a richer 
variety of length scales. 
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ATs = 800 K (Model A1S0) 

ATs = 600 K 

ATs = 400 K 

ATs = 160 K 

ATs = K; Ts = 273 K 

(Model C1S0) 

ATs = K; Ts = 759 K 

(Model H1S0) 
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Fig. 3. — Power spectral density of lateral 
convective flow velocity. Models with differ- 
ent surface temperature contrasts, ATs, an d 
surface temperatures, Ts, are compared. 



5. Rate and Distribution of 
Volcanism 

For surface temperature contrasts 
ATs ^400 K, melt rates show an uneven 
distribution without a coherent global pat- 
tern (Figure Hk, green and orange curves), 
and a total melt production similar to the 
present-day Earth value of ~20 km 3 yr _1 



f lMcKenzie and Bickle I Il988h (Figure gb). 



Above ATs ~400 K, melt occurs preferen- 
tially within upwellings near the hot substel- 
lar point (Figure UK, red and purple curves). 
Due to more persistent day side melting, the 
total melt production rises above the present 
day Earth value by a factor of ~5 for ATs 
= 800 K and ~30 for AT S = 1200 K (Figure 
Hb, red and purple circles, respectively). 





n £ 






> 




XL 
\ 


0.5 


C 


0.4 


o 




3 


0.3 


JO 




istr 


0.2 


Q 






0.1 - 


a) 




m 


0.0 




— ATs=1200K 

— ATs=800 K 
(Model A150) 
ATs=400 K 

— ATs=0 K; 
Ts=273K 
(Model C150) 

— ATs=0 K; 
Ts=759K 
(Model H150) 



-180 -120 -60 60 120 180 

Angular Coordinate [degree] 



0.1 

" 0.05 

O 

E 0.02 

> 

~ 0.01 

a> 

2 0.005 
CO 

H 0.002 
0.001 



: (b) 












: C150 / 


Present Day 




Earth Value 



200 400 600 800 1000 1200 

Surface Temperature Contrast [K] 

Fig. 4. — (a) Distributions of time-averaged 
melt production. Surface temperature con- 
trasts, AT S and surface temperatures, Ts, for 
different models are color-coded as shown in 
the legend. The purple curve peaks at 1.18 m 
kyr -1 . (b) Total melt production, shown el. S 9 
processing rate in mantle volume per billion 
years. The dashed line shows Earth's present- 
day melt production of ~20 km 3 yr _1 , corre- 
sponding to ~0.002 V m Gyr- 1 . 



For planets at the same orbital distance, 
more than double the amount of melt is pro- 
duced in the model with uniform surface tem- 
perature H150 (FigureHb, blue circle) than in 
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the asymmetric model A150 (Figure Hb, red 
circle). Melting is more vigorous for model 
HI 50 because diffuse surface deformation oc- 
curs globally, and is less vigorous for model 
A150 because deformation occurs diffusely 
only near the substellar point. 



6. Discussion and Concluding 
Remarks 

Many recently-discovered exoplanets in- 
habit close-in orbits and this results in high 
(100's to 1000's K) effective temperatures 
(Figure [5|). Planets closer than 0.5 AU 
are likely ti dally locked (Figure EJ black 



dotted line) feasting et al. I Il993l ) and syn- 



chronously rotating around their parent stars. 
At distances of <0.1 AU, Earth-mass plan- 
ets are unlikely to retain an atmosphere, due 
to atmospheric loss eit her by extr eme ultra- 
violet (EUV) heating (|Tian 1120091) or stellar 
wind erosion (jLammer et al. 1 120081 ) (Figure 
purple dashed curves). These boundaries 
move outward for smaller planets that orbit 
around larger stars. Absence of a substan- 
tial atmosphere results in high substellar tem- 
peratures, approximated as T su bst = y/2T e jf, 
and the night side remains cold. Thus, hot 
A-model behavior is expected for close-in or- 
bits, and our results demonstrate that a per- 
sistent hemispheric surface temperature con- 
trast >400 K can maintain a degree 1 pat- 
tern of mantle convection in which the plane- 
tary surface moves preferentially toward sub- 
duction zones on the cold night side. These 
planets should exhibit an inhospitable day 
side with vigorous volcanism and a cold night 



side that allows for more Earth-like tectonics 
and plate-like behavior. Melt production in- 
creases with increasing surface temperature 
contrast, and at 0.13 AU the calculated to- 
tal melting rate is ~5 times higher than for 
a cool planet at 1 AU where plate tectonic 
behavior occurs globally (C-models). 

At intermediate distances (~0.1 AU), an 
atmosphere is more likely to persist and at- 
mospheric redistribution of heat is expected 
to produce a uniform hot surface. For such 
planets, our H-model results demonstrate 
that mobility of the lithosphere prevents the 
formation of coherent plates and the emer- 
gence of Earth-like plate tectonics. Com- 
pared to an asymmetric model at the same 
orbital distance, melt production is enhanced 
by a factor >2, and occurs globally across a 
surface that is characterized by diffuse defor- 
mation. Our results suggest the possbility of 
different feedback mechanisms. Volcanic out- 
gassing is responsible for the formation of 
secondary planetary atmospheres. There- 
fore, asymmetric conditions can only be 
sustained if the atmosphere is continuously 
eroded. If, instead, an atmosphere is re- 
tained, heat redistribution promotes mobile 
lid convection with diffuse deformation and 
further increases global volcanic outgassing. 
This suggests a positive feedback in favor of 
a thick atmosphere, unless another mecha- 
nism modifies the atmospheric balance. For 
example, no melt production is calculated 
for one-plate model planets, a possibility for 
potential close-in equivalents of present-day 
Venus. Atmospheres would be unprotected 
against solar wind erosion if a magnetody- 
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namo is deactivated due to a transforma- 
tion from plate tectonics to a regime that 



core ( 


Buffett 


2002 


2004 




Gaidos et al. 



2010h . 



iner 



Substantial at- £2, 



mospheric loss would allow for a negative 
feedback and such planets could fluctuate be- 
tween symmetric and asymmetric end- mem- 
ber scenarios or reach an intermediate equi- 
librium, depending on the timescale of atmo- 
spheric loss relative to mantle thermal evo- 
lution. The dependence of tectonic regimes 
on o y for our nominal Earth models is con- 
sistent with previous studies performed us- 




spherical geometries (Moresi & Solomatov 


1998 


; Tacklev 


2000al: Richards et al. 


2001 


: Van Heck and Tacklev 2008: 


Folev & Becker 1 


2009J). As in the above- 



Orbital Distance [AU] 

Fig. 5. — Possible occupation of planets in 
orbital distance - stellar mass domain space. 
Purple dotted curves show atmospheric es- 
cape es timates for EUV heating for a 6M E 
planet (ITian 1 120091 ) and stell ar wind erosion 



for a n Earth-sized planet ( iLammer et al. 
20081 ) . Tidal locking is likely at distances 



smaller th an the tidal lock orbi t (black dot- 



ted line, feasting et al. I Il993l )). Colored 



mentioned studies plate-like behavior occurs 
at yield stresses lower than classical esti- 
mates for dry oceanic lithospher^ dots show effective temperatures, T e// , for 
up to ~1000 MPa feohlstedt et al. Ill995l ). Kepler candidate exq planets with R < 2R E 



for this discrepancy 


Moore & Rvmer 


2007; 


Escartin et al. 


2008 


. The weak dependence 



(IBorucki et al. 1 1201 ll ) . For planets with a 



of tectonic regimes on Ra shown in Figure 
[2^ is in agreemen t with numerical results 
of iFoley fc Becker I (120091 ) and s caling laws 
of plat e tectonic convection by iKorenaga 
(j2010al lbl). Because a higher Ra is expected 
for more massive planets, this low sensitivity 
to Ra suggests that our results can be applied 
to planets of various sizes. 



substantial atmosphere, uniform surface tem- 
peratures are expected. Planets without 
an atmosphere can have estimated substellar 
temperatures of T suhsteU ar r ^ j \pXT e jj and a cold 
night side. The theoretical habitable zone lies 
between the 273 and 373 K isotherms where 
liquid water can be maintained. 



Convection in massive "super-Earth" 
mantles may be influenced by their more ex- 
treme pressure and temperature conditions. 
For example, mineral physics calculations 
suggest a viscosity decrease of 2 to 3 or- 
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ders of magn itude for the deep mantles of 
super-Earths (IKarato 11201 lh . in favor of vig- 
orous convective overturn. High mantle pres- 
sure allows for mineral phase transforma- 
tions that do not occur in Earth's man- 



tle (jUmemoto et al. I 120061 ) but which can 
strongly influe nce the dynamics of sup er- 
Earth mantles (IVan den Berg et al. 1120101 ). 



Although important questions remain, 
our simulations demonstrate the strong influ- 
ence that surface temperature contrasts exert 
on mantle convection, surface tectonics, and 
volcanism for close-in rocky exoplanets. Dis- 
tinct scenarios are likely associated with vari- 
ations in albedo, volcanism, and atmospheric 
content that may become astronomically de- 
tectable in the future. 
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